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Abstract 

We define a class of Lorentz invariant Bohmian quantum models for N en- 
tangled but noninter acting Dirac particles. Lorentz invariance is achieved for 
these models through the incorporation of an additional dynamical space-time 
structure provided by a foliation of space-time. These models can be regarded 
as the extension of Bohm's model for N Dirac particles, corresponding to the 
foliation into the equal-time hyperplanes for a distinguished Lorentz frame, to 
more general foliations. As with Bohm's model, there exists for these models 
an equivariant measure on the leaves of the foliation. This makes possible a 
simple statistical analysis of position correlations analogous to the equilibrium 
analysis for (the nonrelativistic) Bohmian mechanics. 

03.65.Bz 
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I. INTRODUCTION 



Among the different approaches to resolving the conceptual problems of quantum theory, 
Bohm's approach is perhaps the simplest. In a nutshell, it consists in adding the most basic 
dynamical variables, obeying additional evolution equations, to the description of a quantum 
system provided by its wave function if>. For nonrelativistic quantum theory the additional 
variables are the positions of the particles, which evolve according to a "guiding equation" 
naturally suggested by the Schrodinger evolution. This theory — usually called Bohmian 
mechanics or the pilot-wave theory — is well understood. It has been analyzed, and its 
connection with the predictions of orthodox quantum theory explained, in the original papers 
of Bohm |l| as well as in later works (see, e.g., One of the main problems remaining 

for the Bohmian (or any other) approach is to find a satisfactory relativistic quantum theory, 
a theory that is fully Lorentz invariant while avoiding the profound conceptual difficulties 
of orthodox quantum theory. 

In his original papers, Bohm had an outline for a "Bohmian" field theory, with fields 
on space-time as the additional variables. A year later he proposed a "Bohmian" model 
for one Dirac particle |5[], which was subsequently extended by Bohm and coworkers to iV 
Dirac particles ||. For this A-particle model the additional variables are, as in Bohmian 
mechanics, the positions k = 1,...,N, of the particles. However, in contrast with 
Bohmian mechanics, the guiding equation for this theory 

dQk _ 4> j cz k ip , y. 

dt ip^ip U 

is ultralocal on configuration space: The right hand side of (|1|) depends only upon the 
value of ip at the positions of the particles and not upon spatial derivatives of ip there. 
Here ip = ^(cjij . . . , q#, t), taking values in the iV-particle spin space (C 4 )®^, solves the 
iV-particle Dirac equation (h = c = 1) 

= S {~ ia k ■ V fc - ecy. k ■ A(q fc , t) 

k=l 

+e$(qjfc,t) + p k m\ip, (2) 

where ct k = {a\, a|, a|), a k = I<g>---®I®a l ®I®---®I, with the z-th Dirac a matrix 
a % at the k-th of the N places, and (3k is defined analogously. $ and A are external elec- 
tromagnetic potentials. (We may of course consider particle-dependent masses m&, charges 
Ck-, and external potentials and A*,, but for simplicity we shall not do so.) We shall call 
this model the Bohm-Dirac model (BD model). Just as with Bohm's proposal for a field 
theory, the BD model requires for its formulation the specification of a distinguished frame 
of reference — in terms of which the actual configuration (Qi, . . . , Qat) and the generic con- 
figuration (qi, . . . , q^v) at time t is defined — and in fact the model is not Lorentz invariant 
if N> 1 g. 

However, for N — 1 this model is Lorentz invariant, and may be formulated in a co- 
variant way: Writing X = X(r) for the space-time point along a trajectory, with (scalar) 
parametrization r, the guiding equation may be written as 
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(3) 



with if; satisfying the Dirac equation 

(ry ■ d — e7 • A — m) ip = 0, (4) 

where 7 • d = and 7 • A = 7 At A M (x). Note that the right hand side of (0), the Dirac 
current j = j M = ip^ip, is the simplest 4-vector that can be constructed from the Dirac 
spinor ip. 

Note also that the parameter r has no intrinsic physical significance, so that equation 
(|) is equivalent to 

dX _ . 
dr 

with arbitrary positive scalar field a = a(x). It is not the field of 4- vectors j (having 
direction and length) that determines the particle motion, but rather the field of directions 
defined by j. In other words, the law for the particle motion could be formulated in a 
purely geometrical manner as the condition that the Dirac current j at every point along 
the trajectory be tangent to the trajectory at that point. 

Because the Dirac current is time-like and divergence free I] 

d-j = 0, 

there is a dynamically distinguished probability distribution on the set of particle paths X{r) 
arising from ([|). Any distribution on this space of paths can be defined by specifying for the 
path the crossing probability for some given equal-time surface So in some Lorentz frame. 
(By this crossing probability we mean the distribution of the point through which the path 
crosses S , which is the same thing as the probability distribution for the position of the 
particle in this frame at the given time.) The distinguished distribution is then defined by 
the crossing probability for So given by g = j° = ip^ip 011 ^0 (with ip suitably normalized), 
which can be written in a covariant manner as j ■ n where n in the future-oriented unit 
normal to the surface. For this distribution the crossing probability for any other equal-time 
surface will also be given by j ■ n, both for the original frame and any other Lorentz frame. 
We may roughly summarize the situation by saying that for the distinguished probability 
distribution, quantum equilibrium holds in all Lorentz frames at all times, with the quantum 
equilibrium distribution given by g = ip^i/j. 

More generally, the crossing probability for any space-like hypersurface S will also be 
given by j ■ n, with n = n(x) the future-oriented unit normal field to S. Moreover, for any 
oriented hypersurface S, the crossing measure (a signed measure that need not be normal- 
ized), which describes the expected number of signed crossings through area elements of S, 



The claims in this and the next paragraph follow directly from the application of the divergence 
theorem (or Stokes' theorem) to an infinitesimally thin tube of paths between So (see below) and 
the relevant hypersurface S. 
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with negatively oriented crossings counted negatively, is, for the distinguished distribution, 
also given by j ■ n, with n = n(x) now the positively oriented unit normal field to S.0 

The iV-particle BD model ([!]) also has a dynamically distinguished probability distribu- 
tion on paths. As a consequence of (fj) g = ip^ip satisfies, in the Lorentz frame in which the 
dynamics is defined, the continuity equation 

! + f> fc -J fc = 0, (5) 

k=l 

where 

Jfc = Qv k = ^at k ijj . (6) 

Thus, if the joint probability distribution for the positions of the N particles is given g = ip^ip 
at some time t = to, then, for the corresponding distribution on paths, it will be given by g = 
at all times t. However, even for this distinguished distribution, quantum equilibrium 
will not in general hold in other Lorentz frames: The joint distribution of crossings of 
equal-time surfaces for other frames will in general not be given by ip^ip' (where ip' is the 
wave function in the relevant Lorentz frame) |7]f|. Nonetheless, Bohm and coworkers have 
argued that the observational content of this model is as Lorentz invariant as the covariant 
formalism of relativistic quantum theory: Since the predictions for results of measurements 
for this model can be regarded as reflected in the configuration of various devices and 
registers — and hence can be derived from probabilities for positions given by g = ip^ip — at 
a common time in the distinguished frame, these predictions must agree with those of the 
usual interpretation. Thus no violation of Lorentz invariance can be detected in experiments 
||. (In particular, the identity of the distinguished Lorentz frame cannot be ascertained by 
means of any possible observation.) 

Lorentz invariance is, however, a delicate issue. Indeed, any theory can be made trivially 
Lorentz invariant (or invariant under any other space-time symmetry), even on the micro- 
scopic level, by the incorporation of suitable additional structure ||. For this reason Bell 
has stressed that one should consider what he has called "serious Lorentz invariance," a 
notion, however, that is extremely difficult to make precise in an adequate way 0. Lacking 
a general criterion, we may nonetheless begin to get a handle on "serious Lorentz invariance" 
by analyzing some specific models. If the models involve additional structure, then whether 
or not we have serious Lorentz invariance will depend, of course, upon the detailed nature 
of this structure. 

In || we have considered a model for which the additional structure for a system of iV 
(noninteracting) Dirac particles is provided by a global synchronization among the particles: 
The trajectories of the particles are such that each one of them at some given space-time 
point is tangent to a vector field determined, given the wave function, by that point and 



2 In this regard it is perhaps worth noting the following: In Minkowski space there is a natural 
duality between divergence-free vector fields and closed 3-forms. Such a vector field defines a 
"deterministic" random path, whose "law" is given directly by the vector field, as in (||), and 
whose statistics are governed by the dual 3-form, in the manner just described. 
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those points along the trajectories of the other particles with which that point has been 
"synchronized." This additional synchronization structure is defined implicitly by the equa- 
tion of motion and the model is not amenable to a statistical analysis in any obvious way. In 
other words, this model is not statistically transparent (see Section IV of ||). Nonetheless, 
even this model provides a counterexample to the widely held belief that a Lorentz invariant 
Bohmian theory for many particles is impossible (unless only product states are allowed). 
In this regard, see also the local model of Squires 0. 

In this paper we shall analyze a statistically transparent counterexample, the "hypersur- 
face Bohm-Dirac model" (HBD model). The basic idea was proposed in (jTO in the context 
of bosonic quantum field theory: In addition to the wave function and field variables, a 
distinguished foliation of space-time — a new element of geometrical structure defining si- 
multaneity surfaces — is suggested as an additional dynamical variable of the theory. These 
surfaces need not be hyperplanes. The defining (Lorentz invariant) equations of the theory 
should describe the evolution of the wave function, the field variables, and the simultane- 
ity surfaces. For a careful philosophical discussion of how this may be compatible with 
some appropriate notion of relativity, even if the simultaneity surfaces should turn out to 
be unobservable, see Maudlin |T] . 



Here we shall consider such a theory, not for fields but for N (noninteracting) Dirac 
particles. We shall discuss an as yet incomplete hypersurface Bohm-Dirac model: The law 
for the evolution of the foliation is not specified, beyond the requirement that it not involve 
the positions of the particles. We present no hypothesis concerning the origin of the foliation, 
but have in mind that the foliation should ultimately be governed by a Lorentz invariant 
law, one that may, for example, involve the iV-particle wave function. (For definiteness we 
shall give some very tentative and less than compelling examples of laws for the foliation in 
Section |CV| .) However, we show in Subsection |111 A| that, regardless of how the foliation is 
determined, the dynamics of the HBD model preserves the quantum equilibrium distribution 
on the leaves of the foliation. Thus the model is amenable to the same sort of statistical 
analysis as for nonrelativistic Bohmian mechanics. This is discussed briefly in Subsection 
|fflB, 



II. THE HYPERSURFACE BOHM-DIRAC MODEL 

A general foliation JF of codimension one on Minkowski space M can approximately be 
thought of as a partition of M into 3-dimensional hypersurfaces. These hypersurfaces are 
the leaves of the foliation. The simplest way to obtain a foliation is by a smooth function 
/ : M — > R without critical points, i.e., df ^ everywhere. The level sets f~ 1 {s) are 
smooth hypersurfaces and form a foliation of M. With the one- form df x , which vanishes 
on the tangent space of the hypersurface through x G M, we may associate by the Lorentz 
metric the normal vector field df(x). If this is time-like everywhere, and thus the foliation 
hypersurfaces space-like, we may normalize df(x) to obtain a unit normal vector field n(x) 
associated with the foliation T . 

We shall consider in this paper only space-like foliations, i.e., foliations by space-like 
hypersurfaces. While obviously different /'s may generate the same foliation J 7 , the future- 
oriented unit normal vector field n is uniquely determined by T . When does a vector field 
v(x) determine a foliation JF such that for all x G M, v(x) is normal to the tangent space 
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of the foliation hypersurface through xl If we denote by V the one-form associated with v 
by the Lorentz metric, then, by Frobenius' theorem, the necessary and sufficient condition 
is that V be completely integrable, V A dV = 0. 

Apart from the foliation, the other dynamical variables of the hypersurface Bohm-Dirac 
model are the usual ones: the wave function if), here for N Dirac particles, and the iV-path, 
the iV-tuple of (everywhere either time-like or light-like) space-time paths, which describes 
the trajectories of the iV Dirac particles. Covariant laws for these dynamical variables 
suggest themselves when we write those of the the Bohm-Dirac model, defined by ([!]) and 
(0), in a coordinate-free, i.e., covariant manner. 

To achieve this we consider first of all the ^-function in the multi-time formalism: For 
iV Dirac particles the wave function ip = ip(x\,X2, ■ ■ ■ ,xn), £ M, takes values in the 
iV-particle spin space (C 4 )®^ and satisfies N Dirac equations 

{ilk ■ d k - e7 fc • A(x k ) - m) ip = , (7) 

k = 1,...,N. Here j k = /<g>---<g>/<g>7<g>/<g>---<g>/, with 7 at the fc-th of the N 
places, and A is an external electromagnetic potential. (Just as with (0), we may of course 
consider particle-dependent masses m k , charges e&, and external potentials A k .) The system 
of equations (|?p is a covariant version of (0); in this multi-time form the Lorentz invariance 
of the law for ip is manifest ||.[] The iV Dirac particles are coupled by the common wave 
function ip. If this is entangled, we have nonlocal correlations between the N particles, 
despite the fact that the particles are noninteracting. 

We shall now develop the guiding law for the iV-path. Note that the numerator of the 
right hand side of ([!]) is given by a current j k , 

j k = ^7?---7fc---7^, 

that involves matrix elements of an operator having as factors the 0-component 7 of a 4- 
vector for all but the k-th particle. Therefore j k can be expressed in a covariant manner by 
replacing 7° in the above expression with 7^ • n, where n is the future-oriented unit normal 
to the t = const hyperplanes, 

jk = V>(Ti • n ) ■ ■ ■ Ik ■ ■ ■ (jN ■ n)i). (8) 

Moreover, the denominator of the right hand side of ([]]) can be expressed covariantly as 
j k ■ n. Then the covariant velocity of the k-th. particle — with respect to the time of a Lorentz 
frame with n as time axis — is 

dXk _ jk 
dt j k ■ n 

Since jk • n = ipdi ■ n) . . . ( r y N ■ n)ip is independent of k, we may reparametrize the paths 
with a parameter s so related to t that t'(s) = jk ■ n to obtain 



3 Note that in the single-time form (^) we can easily add an explicit interaction potential 
V(qi, . . . ,qj\r,t) for the N Dirac particles, while in the multi-time form this is impossible. 
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(9) 



dX k . 

—j— = Jk- 

as 



(10) 



More generally, by further reparametrization, we may obtain dXk/dr = ajk, where a is any 
positive scalar field. The physical particle dynamics — i.e., the N space-time paths defined 
by the equations of motion (and initial conditions) — is invariant under reparametrization. 

A manifestly "parametrization invariant" formulation of the dynamics — that is, such 
that a time parameter plays no role — is easily obtained: The space-time paths for the N 
particles are constrained by the currents jk by requiring that the path for the fc-th particle 
at the point Xk be tangent to the current jk evaluated at Xk and at the the intersection 
points of the paths of the N — 1 other particles with the t = const-hyperplane E 4 containing 
Xk- If we denote by X^(E t ) the intersection point of the path Xk with the hyperplane E t , 
and by XkiT^) a tangent of (or the tangent line to) the path Xk at Xk(T, t ), we may write 
the law for the iV-path as 

X k (E t ) \\j k (X 1 (V t ),...,X N (E t )), (11) 

using the symbol || for "is parallel to." In this geometric formulation the Bohm-Dirac 
dynamics depends upon the Lorentz frame only via its associated foliation into simultaneity 
hypersurfaces £ t , and thus naturally extends to an arbitrary foliation T of Minkowski space- 
time M by curved space-like hyper surf aces :0 

Given such a foliation T and E e J 7 , let X^(E) be the intersection of the path Xk with 
E,0 and let -Xfc(E) be a tangent of (or the tangent line to) the path Xk at X/^E). The law 
of the iV-path X = (X 1; . . . , X^) for the hypersurface Bohm-Dirac model is defined by the 
currents jk naturally extending ([8]) 

jk = i>{li ■ ni) . . . 7 fc . . . (7at ■ n N )?p, (12) 

where rii = n(xi),...,n N = u(xn), with n the future-oriented unit normal vector field 
associated with JF, via the HBD tangency condition (see also Fig. 0) 

X fc (E) || j k {X x {Y),...,X N {Y)). (13) 

[By considering the action of a suitable Lorentz transformation on 7°7 ■ n for arbitrary time- 
like unit vector n (transforming n to (1, 0, 0, 0)), one sees that 7°7 • n is a positive operator 



4 This is in marked contrast with the parametrized dynamics such as given by equations (||) or (|To|). 
which need not extend in anything like the same form to a general foliation since the parametrized 
paths generated by the dynamics need not, in general, respect the foliation. 

5 Note that the paths Xk comprising an iV-path, since they are nowhere space-like, can intersect 
S at most once. This is the main reason why it is important that the foliation T be space-like. 
Of course, also from the physical point of view a synchronization along space-like hypersurfaces 
yields a picture which perhaps makes most sense. We shall assume, without further ado, global 
existence: that a fragment of an iV-path locally satisfying the HBD tangency condition, see (|l3|). 
can be continued in such a manner that each of its paths Xk intersects every S G J-. 
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in spin space C 4 . Hence (7°7i • n\) . . . (7^7^ • n) . . . (7^7^ • tin) is also positive, i.e., j k -n> 
with "=" only if ip = 0. This means that, where it is nonzero, j k is future-oriented and, like 
the path X k , nowhere space-like.] 

We may also write down the equations of motion in the parametrized form analogous to 
(§) or (|IDD. To do so it is convenient to label the hypersurfaces of the foliation using a func- 
tion / : M — > R that generates the foliation as described above, and use this hypersurface 
labeling as the parameter for the particle trajectories — so that Xk(s) is on the hypersurface 
/ -1 (s). From the geometrical characterization of the dynamics fll3|) we know that dX k /ds 
is parallel to j k (Xi(s), . . . , X N (s)), and the scale factor required to ensure f(X k (s)) = s for 
all k and s is easily seen to be l/(df ■ jk)- Therefore 

dX k = j k (Xi(s),...,X N (s)) 

ds df{X k (s))-j k (X 1 (s),...,X N {s)Y 1 } 

For a flat foliation we may choose a Lorentz frame such that the foliation hyperplanes are 
the x° = const-planes, i.e. f(x) = x° for all x. Then n = df = (1, 0, 0, 0) and ( [H] ) reduces 
to the Bohm-Dirac law ([!]). 



III. STATISTICAL ANALYSIS OF THE HBD MODEL 



A. Quantum equilibrium 

We shall show now that for the hypersurface Bohm-Dirac model, with foliation J 7 , there is 
a distinguished probability measure on iV-paths X satisfying the HBD tangency condition 
(fjjD, one for which the distribution of hypersurface crossings Xi(E), . . . , Xjv(S) for £ £ 
T depends only upon ip restricted to E (or, more precisely, to £ ) for ^ satisfying (^). 
We shall say that such a distinguished measure, as well as the corresponding hypersurface 
crossing distribution, is equivariant, defining quantum equilibrium. The physical significance 
of the hypersurfaces £ G T is thus twofold: They serve (via (0)) to define the motion of 
the particles, and, for a quantum equilibrium iV-path, it is "on these hypersurfaces" that, 
manifestly, the "particles are in quantum equilibrium." 

The natural candidate for the equivariant crossing probability density g of the HBD 
model is given by the obvious covariant extension of the equivariant density ip'ij) (= 
tp^ . . . 7^ r '0) °f the BD model: 

g = ^(71 • «i) • • • {lN ■ n N )i). (15) 

To see that this is in fact equivariant, note the following: In view of (|i2|), (i) g = j k - n k and 

j k ■ n k is independent of k. (16) 

Furthermore, (ii) the currents j k are divergence free: 

d k - Jk = 0, (17) 

which follows immediately from (|T2|) using the Dirac equation (0) and its adjoint. These 
two properties of the currents, ( |16"D and (|17D , are the key ingredients for the proof of the 
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equivariance of g. For any current satisfying ( jig ) and (|17|) , for the particle dynamics defined 
by (HP ) Q = jk ■ n k is an equivariant probability density for crossings of the leaves of the 
foliation.^] 

The proof of this assertion consists of two steps: First we determine how an arbitrary 
probability density R on crossings of a foliation hypersurface E evolves under the dynamics 
(|l3l), i.e., we formulate the continuity equation of the hypersurface dynamics. In the second 
step, we show that R = g solves the continuity equation. It then follows that if the proba- 
bility distribution of the "positions of the N particles" on E G T is given by g restricted to 
E, then for any other hypersurface E' G J 7 , the probability distribution of the "positions of 
the N particles" on E' which emerges by transport according to the dynamics ([13]) is given 
by g restricted to E'. Thus g is equivariant. 

Consider thus two infinitesimally close hypersurfaces E and E' belonging to the foliation 
T . The probability distribution of the positions of the N particles on E is given by a density 
i?s : E^ — > R such that 

Prob (particle i crosses E in Sxi, i = 1 . . . N) 
= R s (xi, • • • , x N )5x 1 ■ ■ ■ 5x N . 

By 5x we denote simultaneously an infinitesimal region on E around x and its area (i.e., 
3- volume). Now we compare i?s evaluated at (xi, . . . , x^) G E^ with evaluated at 
(x[, . . . , x' N ) G (S')^, where x' G E' is obtained from x G E via displacement from E to E' 
in the normal direction, see Fig. ^. Let Sx' be the area of the image of the region 6x under 
this correspondence. (Since the projection of the Lorentz metric on E' need not agree with 
the image, under x i— > x', of its projection on E, 5x and 5x' need not agree.) 

Recall from elementary physics that a continuity equation such as (^]) is an expression of a 
local conservation law that, on the infinitesimal level, can be stated as follows: The difference 
between the probability densities i?s on H N and R^ on (E') N (with E' infinitesimally later 
than E) is accounted for by the flux through the lateral sides — to which the hypersurface 
normals are tangent — of the configuration-space-time box between 5x\ x . . . x 5xn C E 
and the corresponding set of (primed) points in (E') , see Fig. U 



N 



Rz>(xi, . . • , x' N )5x[ ■ ■ -8x' N - Rz(xi, . . . , x N ) 8xi ■ ■ -5x N 

.V 

(RsVkjixx, . . . ,x k -i,y,x k +i, • • -,x N ) ■ (u k 5r)(y) dS k , (18) 

t-'fc) 



N r 

8x\ . . . 5xk ■ ■ ■ 5xn j 

fc=l Jd(6x k 



where the ^ on Sx k indicates that this term should be omitted from the product. Here y is 
the integration variable on d(5xk), the (2-dimensional) boundary of 5x k regarded as a region 
in E, dSk is the area element of d(5xk), Uk is the outward unit normal vector field in E to 
d(6xk), Sr(y) is the Minkowski distance between y G E and the corresponding y' G E' (so 
that y' = y + 5r(y)n(y)) and 

v k = ~~~ (19) 
Jk ■ n k 



3 In contrast, the current j k = iplk^ we considered in |8) satisfies (0) but not (|l6| 
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is the covariant velocity of the k-th particle relative to S, see Fig. 



Equation (18) is the continuity equation for the HBD model in the "infinitesimally in- 
tegrated form." It is valid for any hypersurface dynamics defined by ([13]), regardless of 
whether the currents j k satisfy (|l6l) and fll~7|) . However, as we shall now show, if the currents 
do satisfy ( |16| ) and (|iT|), then i?s = f?|s = (jk • H*fc)|s satisfies ([T8|) . 

Sincef] 

g(x[, . . . , x' N )5x[ ■ ■ ■ 5x' N - g(xi, . . . , x N )5x 1 ■ ■ ■ 5x N 
= q(%i, ■ ■ ■ , Xjy)5x^ • • • 5xjy 

Q(^1i * * * 5 X *j\$X\&X 2 ' ' ' &Xjy 

+g(x\, x 2 , . . . , x' N )SxiSx' 2 ■ ■ ■ Sx' N 

-g(x ll x 2 , x' 3 , . . . , x' N )5x 1 5x 2 5x' 3 ■ ■ ■ 8x' N 

H h q(xi, . . . , xjv-;l, x' N )5x! ■ ■ ■ 5x N ^i5x' N 

-g(x 1 ,...,x N )Sx 1 ---Sx N , (20) 
we obtain in this case for the left hand side of flT8"|) (to leading order) 

N 

dxi . . . 8x k ■ ■ ■ 5x N (jk{xi, ...,x' k ,..., x N ) ■ n(x' k )5x' k 



-j fc (xi, . . . ,x k , ■ ■ -,x N ) ■ n(x k )5x k ), (21) 



fc=i 



while the integrand on the right hand side of ( |T8"D becomes (j k -u k )5r dS k . Thus, subtracting 
the right hand side of ( |18|) from ([H]), we obtain (to leading order) the sum over k of the 
integral of j k over the (outward oriented) boundary of the space-time region above 5x k 
between E and £'. But since j k is divergence-free (|T7|), each such term, and hence the sum, 
vanishes. Thus (|T^) is satisfied, establishing the equivariance of g. 

We may also write the continuity equation (|T8"|) in a purely local form: Writing 

6Re(x!, ...,x n ) = Rx>{x[, ...,x' n )- i? s (xi, . . . ,x N ), 

where 

R< s ,(x , 1 ,...,x' N )6xi---5x N 
= Rn'ix^, . . . ,x' N )Sx[ ■ ••8x' N , (22) 



and applying Gauss' theorem to the right hand side of ( jig) 

/ R-zv k ■ u k 5rdS k = div k (Rj:V k 5r(x k )) 5x k 

J dISxi.) 



7 Note that this decomposition is possible because g is defined on M N (with M Minkowski space), 
in contrast with an arbitrary R = (i?s)sG^ r , defined only for iV-tuples belonging to Y< N for some 
S € T , for which therefore such a decomposition is impossible. 



10 



where div^ is the divergence with respect to the k-th coordinate x k on the Riemannian 
manifold £ and is the projection of v k on E, yields 



TV 

5Rz + J2 d ™URzv%5r k ) = 0, (23) 
k=i 

where 5r k = 5r(x k )- 

Using this form we may also check the equivariance of q. To do so, we first "smoothly" 
label the hypersurfaces of the foliation JF by a parameter s € M, increasing in the future 
direction, which may be called a "time parameter," in terms of which (|23|) becomes a stan- 
dard differential equation. The function / : M — > K. that maps any point x G M to the label 
s of the hypersurface E s to which a; belongs generates the foliation in the manner described 
in Section [TT|. In particular, df = \\df\\n, where n is the future-oriented unit normal vector 
field of T . With 5s = \\df k \\5r k , where dfk = df(x k ), we get from (pT9|) that 

Vk5r k = , Jfc ^ = v k Ss, (24) 

with Ufe = dXk/ds the velocity of the fc-th particle in the parametrized formulation of the 
dynamics (III). 

Consider now a coordinate system adapted to our parametrized foliation S s : one coordi- 
nate is clearly given by s, and on one foliation hypersurface we introduce an (arbitrary) coor- 
dinate system p, which is transported to the other foliation hypersurfaces by the flow along 
the normal field, yielding the system of coordinates (s,p), allowing us to write x = (s,p) 
for x e M. Then Xk = (sk,Pk) G S s <^ = s, and the relation between x = (s,p) and 
x' = (s',p f ) from Fig. becomes p = p'. Let Sp be the volume element defined by the 
p-coordinates and let 5x = g(p, s)5p. In these adapted coordinates the continuity equation 
( p3|) assumes, using (|22|) and (H), the more standard form 

^^^W'^O, (25) 



k=l 



with R s (pi, ...,Pn) = -Rs s ((s,pi), . . . , (s,p N )) and g k = g(s,p k ), k = 1, ...,N (and where 
{)^ s is the projection of v k on S s ).0 

For R s = g s , ( pOl) is what lies behind the usual (implication of the) chain rule 

1 d{g x ---g N Q s ) 



91---9N ds 

(26) 



1 d{g{s k ,p k )g{s 1 ,p 1 ,...,s k ,p k )) 



8 This evolution equation depends upon g only through the area-expansion factor arising from the 
normal flow between hypersurfaces, and thus does not really depend upon the choice of coordinates 
on the hypersurfaces. 
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Splitting the 4-divergence into pieces corresponding to variations orthogonal to and varia- 
tions within S s , we obtain 

divj = 110/n ( g f) +div s = (na/irv 53 ') 

where j° is the normal component of j, j° = j ■ n. Setting j = jk and using div jk = ( |TT[ ) 
we then find with ([24]) that 



for all /c. Therefore, in view of (EH), summation over establishes that i? s = £ s satisfies the 
HBD continuity equation (p5|). 



B. Comparison with quantum mechanics 

The statistical analysis of the hypersurface Bohm-Dirac model can be based on the 
assumption that the probability distribution on iV-paths is given by the equivariant density 
g (0) on some simultaneity surface S belonging to the foliation T . Then, by equivariance, 
the statistical predictions of the HBD model (i.e., the crossing probabilities) agree with the 
quantum predictions for positions for any hypersurface in T . But what can be said about 
the statistical predictions concerning a hypersurface which is not part of a member of J-"? 

For one particle the situation is very simple: From the geometrical formulation of the 
HBD model (Section |H[) it follows immediately that the HBD model for one particle is 
foliation-independent, and in fact is the usual one-particle Bohm-Dirac theory given by eqs. 
(||) and (Q), with current j = ipjip. Thus in this case the statistical predictions of the model 
agree with the quantum predictions for position along any hypersurface. 

The situation is analogous for N independent particles: If the wave function ip is a 
product wave function, ip = ^i(xi) • • • i[>n(xn), then it follows from the multi-time Dirac 
equation (0) that ipk satisfies the usual one-particle Dirac equation. Furthermore, the path 
of the k-th particle is tangent to the one-particle current ipkl^Pk and thus independent of 
the paths of the other particles. Moreover g is the product of the corresponding 1-particle 
distributions. Therefore, a product wave function indeed generates a foliation-independent 
motion, the motion of iV independent Bohm-Dirac particles, and we thus have agreement 
with all the quantum position distributions in this case. 

In the general case the situation is more subtle: If the iV-particle wave function is entan- 
gled, it will not in general be the case that the distribution of crossings of hypersurfaces not 
belonging to the foliation agree with the corresponding quantum position distributions f7||| 
(which, in fact, may be incompatible with the crossing statistics for any trajectory model 
whatsoever). However, this disagreement does not entail violations of the quantum predic- 
tions, as has been discussed for the case of the multi-time translation invariant Bohmian 
theory in ||. In fact, insofar as results of measurement are concerned, the predictions of our 
model are the same as those of orthodox quantum theory, for positions or any other quantum 
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observables, regardless of whether or not these observables refer to a common hypersurface 
belonging to 

This is because the outcomes of all quantum measurements can ultimately be reduced to 
the orientations of instrument pointers, counter readings, or the ink distribution of computer 
printouts, if necessary brought forward in time to a common hypersurface in JF, or even 
to a single common location, for which agreement is assured. Nonetheless, this situation 
may seem paradoxical if we forget the non-passive character of measurement in quantum 
mechanics. The point is that for Bohmian quantum theory, measurement can effect even 
distant systems, so that the resulting positions — and hence their subsequently measured 
values — are different from what they would have been had no measurement occurred. 

IV. PERSPECTIVE 

We have presented a hypersurface Bohm-Dirac model for N entangled but noninteracting 
Dirac particles. This model is a covariant extension of the Bohm-Dirac model, which involves 
a foliation by equal-time (flat) hypersurfaces, to arbitrarily shaped (smooth) hypersurfaces. 
How natural is this model? 

When looking for a relativistic extension of nonrelativistic Bohmian mechanics one in- 
evitably encounters two central, very different problems: that such an extension must involve 
a mechanism for nonlocal interactions between the particles, and that quantum equilibrium 
cannot hold in all Lorentz frames. For both of these problems the additional space-time 
structure provided by a foliation yields the most obvious solution: The motion of each par- 
ticle at a point x e M depends upon the paths of the other particles via the points at which 
they intersect the leaf of the foliation containing x, and we have an equivariant density on 
the leaves of the foliation. 

And the simplest way to achieve this, in a covariant manner, for a Dirac wave function 
if), is via the current fll2f) : Form the natural tensor ipji . . . evaluated at x and the other 
intersection points, and contract in the slots corresponding to the other particles with the 
N — l unit normals to the hypersurface at the corresponding points, to obtain the divergence- 
free 4-vector jf., the tangent to the trajectory at x. Thus, the dynamics of the HBD model is 
the simplest Lorentz invariant dynamics compatible with the structure at hand, namely, the 
Dirac wave function and the foliation. Furthermore, the simultaneous normal component 
Q = jk ' n k is an equivariant density on the leaves of the foliation. 

It should be stressed, however, that the Lorentz invariance of the HBD model is — in Bell's 
sense — "serious" only if the foliation is regarded as an additional objective dynamical — in 
contrast to absolute — structure in the theory (and in the world, if the theory is to describe 
the world). It is this structure that is the innovation of what has been proposed here and 



9 This conclusion requires the rather dubious assumption that the relevant measurements can be 
understood in terms of noninteracting Dirac particles. However, in order to talk coherently about 
the quantum predictions for a model, it must be possible to understand measurement processes 
in terms of that model. The remarks we are making here would also be appropriate for the more 
realistic models for which this would be true. 
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Ill 



T0|1 , n °t the model per se, which is indeed a rather straightforward covariant extension 
of the BD model. 

However, in this paper we shall not try to find a "serious" law for the foliation T or, what 
amount to the same thing, its normal vector field n. As a toy example, however, the foliation 
law could be given by an autonomous equation for n, such as d v n^ = 0. Another class of 
toy examples involves a vector field n constructed from the wave function ip(xi, . . . ,xn)'- 
Consider the space-time vector fields v^(x) = (^7^)(fl, . . . , xfj^i, x, x~i+i, . . . , x~n), where 
(xi, . . . , x/v) is a point fixed in a Lorentz invariant way, for example as a maximum of ipip. 
(Simply considering v%(x) = (ip^ip) (x , . . . , x) is not a good idea, since this will be zero for 
antisymmetric (fermion) wave functions.) Now one may set n equal to the integrableQ part 
of some vm- 

A further possibility, which may be more serious, is to have, in addition to the particle 
degrees of freedom, an independent quantum field M that determines the foliation. Assume 
that for any quantum state $ of the field, ($, <fi^) is time-like and completely integrable. 
Then for any state \I/ of the particle-field system, set n M = (ty, (p^). Suppose that the 
particle and the field degrees of freedom are both dynamically and statistically independent, 
i.e., that there is neither quantum interaction nor entanglement between these degrees of 
freedom, so that in particular the full wave function \I/ = ip ® Then we may define the 
foliation by the normal field n M . The 0-field can be regarded as very roughly analogous to 
a Higgs field, producing a kind of spontaneous symmetry breaking, where by choice of $ a 
particular foliation is determined, and relativistic invariance thereby broken. 
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10 For an arbitrary vector field v ll (x), the Fourier transformed v^(k) may be split into v^(k) = 
v u {k)k v k^ / {k\k x ) and v^_(k) = v^{k) — v^(k). The inverse Fourier transformed v^(x) satisfies the 
integr ability condition d^v^ u — d u v^ = 0. 
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FIGURES 




FIG. 1. Geometrical formulation of the dynamics for a system of three particles: For each 
particle the path of that particle, say particle 1 at xi, must be tangent to the 4- vector j\ 
which is determined by: 1) the intersections X2 and X3 of the trajectories of the other two 
particles with the hypersurface S containing xi, 2) the future-oriented unit normals n<2 and 
723 at these points, and 3) the wave function of the system evaluated at x\, X2 and £3: 
ji = i/>(xi,X2, £3)71(72 • ^2X73 • n 3 )ip(xi,x 2 ,x 3 ). 
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FIG. 2. Motion of two particles in one space dimension from hypersurface £ to £': space-time 
view. We have indicated the positions of the primed points x' k obtained from Xk via displacement 
from £ to £' in the normal direction, and the images 5x' k of the regions 5xt under this correspon- 
dence. The point on £' to which particle k moves when starting at x^ G £ is given (to leading 
order) by x& + v^St^ with = jk/(jk • n k), where St^ is the Minkowski distance between x^ and 



FIG. 3. Conservation of probability for a system of two particles in one space dimen- 
sion: configuration-space-time view with, for simplicity, the hypersurfaces drawn straightened 



boundary points of Sx^.) The change of the probability of particle 1 being in 5x\ and par- 
ticle 2 being in 8x2 from hypersurface £ to £' is accounted for by the single particle fluxes 
through the lateral sides of the configuration space time box between 5x\ x 6x2 C £ 2 and 
the corresponding set of primed points on (£') 2 i i.e., R^'(x[, x' 2 ) Sx'-^Sx^ — Ry,(xi,X2)5x\5x2 = 

- ((RxVi)(yi,X2) • 0i5r)(yi) + (R^v 1 )(y 1 ,x 2 ) ■ («i5r)(yi)) Sx 2 

- ((RT,v 2 )(xi,y2) ■ (u2ST)(y 2 ) + (-Re^K^i, 2/2) • («2^t)(2/2)) Sx\. Eq. (jl|) is the natural extension 
of this formula to N particles in Minkowski space. 





8x'i 
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